arXiv:hep-ph/0211303v2 20 Nov 2002 


The Resummed Photon Spectrum in Radiative T Decays 
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We present a theoretical prediction for the photon spectrum in radiative T decay including the 
effects of resumming the endpoint region, —> Mr/2. Our approach is based on NRQCD and the 
soft collinear effective theory. We find that our results give much better agreement with data than 
the leading order NRQCD prediction. 


The radiative decay T ^ Xj was first investigated 
about a quarter century ago Q. The conventional wis¬ 
dom at was that this process is computable in perturba¬ 
tive QCD due to the large mass of the b quarks. Since 
then, we have learned much about quarkonium in gern- 
eral and this process in particular [| S 0 . In 
addition, CLEO is currently taking data on the low ly¬ 
ing T resonances and will soon be able to update their 
original measurement of this decay |^. It is thus timely 
to reexamine the theoretical predictions for this rate. 

The current method for calculating the direct radiative 
decay of the T is by using the operator product expan¬ 
sion (OPE), with the operators scaling as some power of 
the relative velocity of the heavy quarks, v, given by the 
power counting of Non-Relativistic QCD (NRQCD) [^. 
The u ^ 0 limit of NRQCD coincides with the color- 
singlet (CS) model calculation of [^. 

This picture is only valid in the intermediate range of 
the photon energy (0.3 z < 0.7, where z = 2Ej/M, 
and M = 2mb). In the lower range, z < 0.3, photon- 
fragmentation contributions are important i I- At 
large values of the photon energy, z > 0.7, both the per¬ 
turbative expansion Q] and the OPE Q break down. 

The breakdown at large z is due to NRQCD not includ¬ 
ing collinear degrees of freedom. The correct effective 
field theory is a combination of NRQCD for the heavy 
degrees of freedom and the soft-collinear effective theory 
(SCET) 1^, |W| for the light degrees of freedom. In a pre¬ 
vious paper^l we applied SCET to the color-octet (CO) 
contributions to radiative T decay. Here we treat the CS 
contribution at the endpoint within SCET. In this letter 
we present the main results of the analysis, and leave the 
details to a companion paper ||l^ . 

The inclusive photon spectrum can be written as a sum 
of a direct and a fragmentation contribution 0 , 
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where in the direct term the photon is produced in the 
hard scattering, and in the fragmentation term the pho¬ 
ton fragments from a parton produced in the initial hard 
scattering. The fragmentation contribution has been well 
studied in Ref. Q] , and we do not add anything new here. 

The direct contribution can be calculated using the 


OPE, where the rate can be written as 

^ = ^C„(M,z)(T|a|T). (2) 
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The Ci are short-distance coefficients, calculable as a per¬ 
turbative series in as{M), and the O are NRQCD oper¬ 
ators, scaling with specific powers of u <C I. 

At leading order in v only a CS term contributes. 
The CS operator 0(1, ^^i) creates and annihilates a CS 
quark-antiquark pair in a ^Si configuration, and is mul¬ 
tiplied by the CS coefficient, which at leading order is 
proportional to a‘^{M). There are also CO contributions 
down by v'^. Two of these, proportional to the CO 
and ^Pq matrix elements (MEs), give rise to large en¬ 
hancement at the endpoint [Q. Since the CO ^5*0 and 
^Pq MEs are unknown, and since the data does not show 
any enhancement near the upper endpoint, we set them 
to zero. However, we include the CO ME. It has a siz¬ 
able fragmentation contribution, but becomes negligible 
as z increases, and thus does not conflict with data [||. 

The lowest order CS direct rate is Q 
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and Cfe = —1/3. The ME is related to the wavefunction 


(T|Oi(35i)|T) = ^|i?(0)r 
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The Os correction to this rate was calculated numerically 
in Ref. 0 , leading to small corrections over most of phase 
space. In the endpoint region, however, the corrections 
are of order the leading contribution. 

In the endpoint region, the outgoing gluons are moving 
back-to-back to the photon, with large energy and small 
invariant mass (ie, a collinear jet). We must, therefore, 
couple NRQCD to SCET [^. The scales, set by the 
lightcone momentum components of the collinear parti¬ 
cles, are widely separated. If we choose p~ to be 0{M), 












then pj_/p ~ A, and /p ~ A^, where A is a small 
parameter. Here the collinear scale is 
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Pc M\/l — z ^ ^ MAqcd . (6) 

Thus A is of order — z ~ sfkqpajM. 

There are two types of fundamental objects in SCET 
(fields and Wilson lines) and two separate sectors 
(collinear and usoft). In the collinear sector there is a 
fermion field a gluon field g, and a Wilson line 
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Collinear fields are labeled by a direction and the large 
components {n ■ q, q±). The operator projects out the 
momentum label. Likewise in the usoft sector there is a 
fermion field g„s, a gluon field A!^^, and a Wilson line Y. 
Operators are constructed out of these objects such that 
they are gauge invariant. Thus, operators with collinear 
gluons are built out of the homogeneous (order A) com¬ 
ponent of the collinear field strength, PB^ = [ p^ , 



FIG. 1: Feynman diagram for the leading order jet function. 
Collinear gluons are represented by a spring with a line. 


where the forward scattering amplitude T{z) is 

T{z) = -ij <Pxe-^‘i-^T\TJl{x)UQ)\T)g^^ . (13) 

The T indicates time ordering. Matching onto SCET the 
forward scattering amplitude can be written as 

T{z) = ^H{uj,p)Tes{uj,z,p), (14) 

UJ 

where 

Tesii 0 ,z,p)= J d£+Me+ +Mil- z)]Sie+). (15) 


B^^ = ^W\V^+gs{A>:,^g)±)W. ( 8 ) 
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We now write down the leading operator. Aside from 
B±, we also need the NRQCD heavy quark and antiquark 
fields, tjjp and X-p, which transform only under usoft 
(not collinear) gauge transformations. A CS ^Si bb pair 
decays into a photon and a colorless jet of gluons. We 
must, therefore, include two of the B± fields in a colorless 
configuration, and the only operator is 


0(1,3^1) = (9) 

xLp^VpTr{H2 iPt) , 

where P^ acts to the left. Momentum conservation forces 
the momentum of the jet to be M, so B'J_{P + P^)B^ = 
—MB^B^. Introducing V- = P — P\ Eq. (^) becomes 

0{lpSi){M) = (10) 

xlpCT^PpYr{BlT%l^\M,V.)Bi] . 

Matching onto QCD at tree level, we obtain 
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for a transverse photon, where n ■ q- = n ■ q — n ■ q' and 
= gm' _ + rYn^)l2. 

The inclusive T ^ Xy rate can be factored into hard, 
jet, and usoft functions at the endpoint. Using the optical 
theorem the inclusive spectrum can be written as 
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After decoupling usoft degrees of freedom [^, the CS jet 
function is defined as 

(0|TTr[H®“<5„.P_H®^] (x)Tr[H®“'^,,-,p_(0)|0) 

+9f 9^)5.,.' I ^e-*'=-J.(fc+), 

(16) 

and the CS usoft function is defined as 

5(^+) = J (17) 

X {T\T[)Pl(TiX-p] (a;") [xlp-cTi'i/'p'] (0)|T) 

= (T|V'p<7'»X-p^(w • d - £~^)x-p''^^'^P'W ■ (18) 

The hard coefficient H(ijj,pl) can be calculated perturba- 
tively in an expansion in as{M). At tree level we obtain 
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At the collinear scale pc we perform an OPE, integrate 
out collinear modes and match onto a non-local usoft 
operator, Eq. (|l7|), convoluted with a Wilson coefficient, 


T{z) = J di+S{i+)nj[£+ + M(1 - z)]. (20) 

To leading order in as{My/l — z), the jet function is 
calculated from the Feynman diagram shown in Fig. |^. 
Evaluating the diagram gives 
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FIG. 2: Diagrams needed to calculate the CS counterterm. 
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[(fc+/M)(i-e)]^ 

and taking the imaginary part we obtain 

r*l 




ImJ^(A:+) = ^0(fc+)y d^Su;,Mi- (22) 


Combining, we get 
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where the 2 M/M^ accounts for the non-relativistic nor¬ 
malization of the T state in the usoft function. This is 
precisely in form given in Eq. (|^), and it is straightfor¬ 
ward to read off Hj. 

Using Eq. (|l^, we can integrate over £+, giving 
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x(T|V'^criX_p0[w -8 + M{1 - 2 :)]xLp/Cri^/)p.|T) 
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where we used the results of Ref. Q for the final line. 
Plugging into Eq. ( p^ gives the z —> 1 limit of Eq. (i- 
At this point, large logarithms will appear in the jet 
function at higher order. This can be avoided by running 
operators from M to /ic, which sums logs of 1 —z. To run 
the CS operator, we calculate the counter term, deter¬ 
mine the anomalous dimension, and use this in the reno- 
malization group equations (RGEs). The graphs needed 
are shown in Fig. Diagrams involving usoft gluons van¬ 
ish. Feynman rules for the vertex operators are given in 


Ref. 1^. We perform our calculation in Feynman gauge, 
and obtain a relatively simple result for the one-loop UV- 
divergent term 
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This depends on the large momentum component of the 
gluons. The divergent piece must be canceled by the 
counterterm Za/Zo — 1 , where Zq is the CS vertex coun¬ 
terterm, and Z 3 is the gluon wavefunction counterterm 


Z, = 1 
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The anomalous dimension is obtained through the 
standard method, and the RGE for the coefficient is 

^ 7 (^,a;)r(i-'^i)(^,a;). (28) 

d/i 


Solving this equation gives 
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Logarithms of the form hi{fj./M) have been summed into 
r(i. and any logarithms in the operator are of 

the form ln(/ic//r.). If we take /r ~ /ic all large logarithms 
of the ratio /ic/M will sit in the coefficient. 

We now obtain the resummed rate, by substituting 
Eq. (H) into Eq. (|2|), giving 
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Again integrating over £'^ and inserting into Eq. 
resummed CS contribution to the decay rate is, 
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FIG. 3: The inclusive photon spectrum, compared with data 
The theory predictions are described in the text. 

We can expand in as{M) to obtain an analytic expression 
for the next-to-leading logarithmic contribution 

= zi^l + ^[CA{2'K^-17) + 2nf] ln(l-z)|. (33) 

As z approaches one the 0{as) term becomes of order 
one, precisely the behavior observed in Ref. |Q|. The 
resummed result does not suffer from this problem. 

We now combine the different contributions to obtain 
a prediction for the photon spectrum. We will marry our 
expression for the CS spectrum in the endpoint with the 
leading order result by interpolating between the two 

1 dTint / 1 \ 1 dr,esum , , 

To dz VTo dz To dz ' ^ ’ 

Before we proceed we need the NRQCD MEs. We can 
extract the CS ME from the T leptonic width. The CO 
MEs are more difficult to determine. NRQCD predicts 
that the CO MEs scale as compared to the CS ME. 
In Ref. [|^ it was argued that an extra factor of l/2Nc 
should be included. We set the and ^Pq MEs to zero, 
and the ^Si ME to (T|C>8(^5'i)|T) = z;4(T|C>i(3S'i)|T), 
where we use = 0.08. 

The CLEO collaboration measured the number of pho¬ 
tons in inclusive T(IS') radiative decays ||^. The data 
does not remove the efficiency or energy resolution and 
is the number of photons in the fiducial region, | cos0| < 
0.7. In order to compare our theoretical prediction to the 
data, we integrate over the barrel region and convolute 
with the efficiency that was modeled in the CLEO paper. 
We do not do a bin-to-bin smearing of our prediction. 


In Fig. ^ we compare our prediction to the data. 
The error bars on the data are statistical only. The 
dashed line is the direct tree-level plus fragmentation re¬ 
sult, while the solid curve includes the resummation in 
Eq. (|34|). For these two curves we use the Us extracted 
from these data, as{Mr) — 0.163, which corresponds to 
as{Mz) = 0.110 1^. The shape of the resummed result is 
much closer to the data than the tree-level curve, though 
it is not a perfect fit. We also show the Eq. (^) plus 
fragmentation result, using the PDG value of as{Mz), in¬ 
cluding theoretical uncertainties, denoted by the shaded 
region. To obtain the darker band, we first varied the 
choice of mb between 4.7 GeV < < 4.9 GeV and 

the value of as within the errors given in the PDG, 
as{Mz) = 0.1172(20) |Q. Varying mb and as mod¬ 
ifies the extraction of the GS ME from 3.31 GeV^ to 
3.56 GeV^. We also varied the collinear scale, pic from 
— z)/2 < pLc < Myj2{l — z). Finally, the lighter 
band also includes the variation, within the errors, of the 
parameters for the quark to photon fragmentation func¬ 
tion extracted by ALEPH 1^. The low z prediction is 
dominated by the quark to photon fragmentation coming 
from the GO channel. We did not assign any error 
to the CO ME. Since it is unknown, there is a very 
large uncertainty in the lower part of the prediction that 
we decided not to show. Note that the GO ^S'o and ^Pq 
contribution increases the theoretical prediction at the 
upper endpoint ||^. It is thus clear the data favors a 
very small value for the CO ^S'o and ^Pq MEs. This is 
why we set these to zero in our analysis. Negative values 
for these MEs are possible, and would give a bit better 
fit to the shape. 
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